High-accuracy critical exponents of 0{N) hierarchical sigma models 



in 
o 
o 

(N 

> 
O 



> 



Ap. 



J. J. Godina 
Dep. de Fisica 
CINVESTAV-IPN 
Post. 14-740, Mexico, D.F. 



07000 



L. Li and Y. MeuriccQ 

Department of Physics and Astronomy 
The University of Iowa 
Iowa City, Iowa 52242, USA 

M. B. Oktay 
School of Mathematics 
Trinity College 
Dublin, Ireland 
(Dated: February 1, 2008) 

We perform high-accuracy calculations of the critical exponent 7 and its subleading exponent for 
the 3D 0{N) Dyson's hierarchical model for A'^ up to 20. We calculate the critical temperatures 
for the nonlinear sigma model measure 5{(f).4> — 1). We discuss the possibility of extracting the first 
coefficients of the 1 /N expansion from our numerical data. We show that the leading and subleading 
exponents agree with Polchinski equation and the equivalent Litim equation, in the local potential 
approximation, with at least 4 significant digits. 
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The large N limit and the 1/N expansion Q, |^ Q 
appear prominently in recent developments in particle 
physics, condensed matter and string theory 0,|SIEI3- 
For sigma models, the basic gap equation can be obtained 
by using the method of steepest descent for the functional 
integral 0,0- For N large and negative, the maxima of 
the action dominate instead of the minima and the radius 
of convergence of the expansion should be zero. In 
order to turn a expansion into a quantitative tool, 
we need to: 1) understand the large order behavior of 
the series, 2) locate the singularities of the Borel trans- 
form and, 3) compare the accuracy of various procedures 
with numerical results for given values of N . Calculat- 
ing the series or obtaining accurate numerical results at 
fixed N are difficult tasks and we do not know any model 
where this program has been completed. For instance for 
the critical exponents in three dimensions, we are only 
aware of calculation up to order 1 /N'^ in Ref . [3, flOl fllj . 
Several results related to the possibility (or impossibil- 
ity of resumming particular 1/N expansions are known 
El El El. Overall, it seems that there is a rather pes- 
simistic impression regarding the possibility of using the 
expansion for low values of N. For this reason, it 
would be interesting to discuss the three questions enu- 
merated above for a model where we have good chances 
to obtain definite answers. Dyson's hierarchical model 
0,01 is a good candidate for this purpose. 



In this Brief Report, we provide high- accuracy numer- 



ical values for the critical exponent 7, the subleading 
exponent A and the critical parameter [3c for the 3Z3 
0{N) hierarchical nonlinear sigma models. These quan- 
tities appear in the magnetic susceptibility near Pc in the 
symmetric phase as 

X=(/3e-/5)-^(Ao-H^i(/3c-/3)^ + ...) • (1) 

The method of calculation of the critical exponents 
used here is an extension of one of the methods described 
at length in the case of iV = 1 E3 ^-nd will only be 
sketched briefly. On the other hand, the accuracy of 
the approximations used depend non trivially on N as 
we shall discuss later. The RG transformation can be 
constructed as a blockspin transformation followed by a 
rescaling of the field. For Dyson's hierarchical model, the 
block spin transformation affects only the local measure. 
The RG transformation can be expressed conveniently in 
terms of the Fourier transform (denoted R hereafter) of 
this local measure. In the following, we keep the 0{N) 
symmetry unbroken and the Fourier transform will de- 
pend only on k.k = u. Here fc is a source conjugated 
to the local field variable 0. Replacing fc by w and the 
second derivative by the TV-dimensional Laplacian in Eq. 
(2.5) of Ref. [13, we obtain the RG transformation for 
the Fourier transform of the local measure: 
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K.^(.)oce[-^K^«^+^^^)] iR,,^icu/4)f , (2) 

where c = 2^^^/-^ in order to reproduce the scaling of a 
Gaussian massless field in D dimensions. D = 3 here- 
after. We fix the normalization constant by imposing 
i?„jv(0) = 1 so that i?,i.Ar(fc) has a simple probabilistic 



2 



interpretation 17]. In the following, the calculations will 
be performed using polynomial approximations of degree 



Rn,N{k) ~ 1 + a„4ii + a„,2U^ H h 



(3) 



The finite volume susceptibility for 2" sites is related to 
the first coefficient by the relation Xn = — 2a„^i(2/c)". 
The truncated recursion formula for the a„ m reads 



(4) 



with 



_ r(^ + i)r(z + 7V/2) 1 ^c^,^ 

" r(m + l)r(m + iV/2) (/-m)!M^ ^ '^^ 

(5) 

We emphasize that in the above formula and in our nu- 
merical calculations, no truncation is applied after squar- 
ing and so the sum in Eq. I^J does extend up to 2lmax- 
Since the derivatives appear to arbitrarily large order in 
Eq. and can lower the degree of a polynomial of order 
larger than Imax, this affects all the coefficients of order 
less than Imax- This procedure has been discussed and 
justified in Ref. [l^ . 

The critical exponents appearing in Eq. are ob- 
tained by calculating the eigenvalues Ai,A2,... of the 
matrix da„+i^i/dan_rn at the nontrivial fixed point. The 
exponents 7 and A, can be expressed as 



7 = 



ln(2/c) 
In(Ai) 



A 



ln(A2) 



In(Ai) 



(6) 



The critical exponents are universal and, within nu- 
merical errors, independent of the manner that we ap- 
proach the nontrivial fixed point. In the following, we 
have mostly started with the local measure of the non- 
linear sigma model d{(f>.(j)~l). The corresponding Fourier 
transform reads 



(-i)'^'r(^) 
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A motivation for this choice is that, as we will explain 
below, the value of Pc can be calculated in the large N 
limit. Other measures have also been used in order to 
check the universal values of the two exponents. 

The asymptotic behavior of the ratio a„+i^i/a„.i al- 
lows us to decide unambiguously if we are in the sym- 
metric phase (where the ratio approaches c/2 ~ 0.63) 
or in the broken phase (where the ratio approaches c). 
Using a binary search, one can determine the critical 
value of (3 with great accuracy. As this critical value de- 
pends on Imax, we denote it Pdlmax)- When Imax 00, 
(3c{lmax) Pc- The rate at which this limit is reached 
depends on N. This is illustrated in Fig. ^ where we 
see that in order to reach (3c with a given accuracy, we 



need to increase Imax when N increases. In Fig. 12 we 
give the minimum Imax necessary for Pdlmax) to share 
20 significant digits with f3c- Imax — 22-1- 6.2iV°-^ is a 
good fit for Fig. 
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FIG. 1: logi 
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calculated for /max ~ 40 to In 



= 10 

60 for iV = 10 (filled circles), iV = 11 (empty circles), N = 12 
(empty triangles) up to TV = 30 (empty squares). 

The nontrivial fixed point for a given value of Imax 
can be constructed by iterating sufficiently many times 
the RG map at values sufficiently close to /?c(^max)- In 
order to get an accuracy e for the fixed point for that 
value of Imax , we need to iterate n times the map until 



A5 



(8) 



in order to get rid of the irrelevant directions. At the 
same time, we want the growth in the relevant direction 
to be limited, in other words. 



|/3-/?c(/max)|A? <e 



(9) 



Combining these two requirements together with Eq. © 
we obtain 



|/?-/3c(;max)| 



(10) 
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FIG. 2: Minimal value of Zmax in order to have 

l/3c(imax)-/3c(oo))| on AT 

' " ' — ^ = — 20 versus N . 
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TABLE I: (3c and the first two eigenvalues for A'' = 1 ... 20. 



iV 




Ai 




1 


1.1790301704462697325 


1 


427172478 





8594116492 


2 


2.4735265752919854000 


1 


385743490 





8563409066 


3 


3.8273820333573397671 


1 


354668326 





8506945150 


4 


5.2111615635533656165 


1 


332749866 





8440522956 


5 


6.6104153462855068435 


1 


317578283 





8376436747 


6 


8.0181114053706725941 


1 


306955396 





8320345022 


7 


9.4307096447427796882 


1 


299321025 





8273378172 


8 


10.846330737925124699 


1 


293666393 





8234676785 


9 


12.263918029354988652 


1 


289354227 





8202833449 


10 


13.682844072802585664 


1 


285978489 





8176485461 


11 


15.102717572108367579 


1 


283274741 





8154492652 


12 


16.523283812777939366 


1 


281066141 





8135953137 


13 


17.944370719047342283 


1 


279231192 





8120168555 


14 


19.365858255947423937 


1 


277684252 





8106600963 


15 


20.787660334686062513 


1 


276363511 





8094834857 


16 


22.209713705054412233 


1 


275223389 





8084547150 


17 


23.631970906283518487 


1 


274229622 





8075484440 


18 


25.054395659078177206 


1 


273356000 





8067446107 


19 


26.476959772907788848 


1 


272582158 





8060271793 


20 


27.899641020779716433 


1 


271892050 





8053832116 



This is an order magnitude estimate, however it works 
well except for iV=l where we need to pick j3 slightly 
closer to the critical value. By "working well" , we mean 
that if we go closer to the critical value, changes smaller 
than e are observed in the first two eigenvalues. The 
numerical results for e = 10"^" and N up to 20, are 
given in the Tables HI and ITll for the values of Imax of Fig. 
121 Errors of 1 or less in the last printed digit should be 
understood in all the tables. 

As N increases, the values displayed in Table HH seem 
to slowly approach asymptotic values. This is expected. 
Using the general formulation of Ref. yja together with 
the particular form of the propagator [l9j for the model 
considered here, one finds the leading terms 

7 ~ 2 + ai/iV+... (11) 
A ~ \ + hi/N + ... 
Pc/N ~ (2-c)/(2(c-l))+ci/iV + ... . (12) 

The magnitude of the coefficients oi, 6i, ci of the lead- 
ing 1/iV corrections can be estimated by subtracting the 
asymptotic value and multiplying by N . The results are 
shown in Table UTTl They indicate that oi ~ —1.6, hi ~ 
—2.0, ci ~ —0.57. It seems possible to improve the accu- 
racy by estimating the next to leading order corrections 
and so on. However, the stability of this procedure is 
more delicate and remains to be studied with simpler 
examples. 

We now compare the exponents calculated here with 
those calculated with three other RG transformations 
[20I I21I [2^ . As we proceed to explain, the exponents 



TABLE IL 7, A and P^/N for A = 1 . . . 20. 



N 


7 


A 




1 


1.29914073 


0.425946859 


1 


179030170 


2 


1.41644996 


0.475380831 


1 


236763288 


3 


1.52227970 


0.532691965 


1 


275794011 


4 


1.60872817 


0.590232008 


1 


302790391 


5 


1.67551051 


0.642369187 


1 


322083069 


6 


1.72617703 


0.686892637 


1 


336351901 


7 


1.76479863 


0.723880426 


1 


347244235 


8 


1.79469274 


0.754352622 


1 


355791342 


9 


1.81827105 


0.779508505 


1 


362657559 


10 


1.83722291 


0.800424484 


1 


368284407 


11 


1.85272636 


0.817977695 


1 


372974325 


12 


1.86561092 


0.832855522 


1 


376940318 


13 


1.87646998 


0.845589221 


1 


380336209 


14 


1.88573562 


0.856588705 


1 


383275590 


15 


1.89372812 


0.866171682 


1 


385844022 


16 


1.90068903 


0.874586271 


1 


388107107 


17 


1.90680338 


0.882027998 


1 


390115936 


18 


1.91221507 


0.888652409 


1 


391910870 


19 


1.91703752 


0.894584429 


1 


393524199 


20 


1.92136121 


0.899925325 


1 


394982051 


00 


2 


1 


2(c- 


jy = 1.42366.. 



should be the same in the four cases (including ours). 
The change of coordinates that relates the RG trans- 
formation considered here and the one studied in Ref. 
|22J is given in the introduction of [23 (for L ^ 2^/^). 
The fact that the limit L — > 1 in the formulation of Ref. 
[2^ yields the Polchinski equation in the local potential 
approximation studied in Ref. [IJ is explained in Ref. 
|24l | . Consequently, these two RG transformations should 
be the same in the linear approximation. Finally, Litim 
[23, 113 proposed an optimized version of the exact RG 
transformation and suggested [2^ that it was equivalent 
to the Polchinski equation in the local potential approx- 
imation. The equivalence was subsequently proved by 
Morris j27i |. 

To facilitate the comparison, we display = 7/2 (since 
77 = here) and uj = A/v in Table Hvl Our results 
coincide with the 4 digits given in column (2) of Table 3 
(for v) and 4 (for lo) in |23|. They coincide with the six 
digits for v given in the line d = 3 of Table 8 of |23| for 
N= 1, 2, 3, 5 and 10. However, we found discrepancies 
of order 1 in the fifth digit of v and slig htly larger for lu 
with the values found in Table 1 of 2U\. Our estimated 
errors are of order 1 in the 9-th digit. For N = 1, this is 
confirmed by an independent method For N = 2, 

3, 5, and 10, this is confirmed up to the sixth digit [2^ . 
Consequently, a discrepancy in the 5-th digit cannot be 
explained by our numerical errors. Note also that for 
N > 2, a is more negative than for nearest neighbor 
models [TH . 

In summary, we have provided high-accuracy data for 
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7, A and (3c for N up to 20. It seems likely that a few 
terms of the 1/N expansion for these three quantities can 
TABLE III: N{2 - 7), iV(l - A) and N{ ^^~_\^ - ^) for iV = be estimated from this data. Work is in progress to cal- 
1 ... 20.. culate these expansions independently by semi-analytical 

methods and learn about the asymptotic behavior of the 



iV 


N{2 ~ 7) 


N{1 - A) 


V2(c-1) 


1 


0.7009 


0.5741 


0.2446 


2 


1.167 


1.049 


0.3738 


3 


1.433 


1.402 


0.4436 


4 


1.565 


1.639 


0.4835 


5 


1.622 


1.788 


0.5079 


6 


1.643 


1.879 


0.5239 


7 


1.646 


1.933 


0.5349 


8 


1.642 


1.965 


0.5430 


9 


1.636 


1.984 


0.5490 


10 


1.628 


1.996 


0.5538 


11 


1.620 


2.002 


0.5576 


12 


1.613 


2.006 


0.5606 


13 


1.606 


2.007 


0.5632 


14 


1.600 


2.008 


0.5654 


15 


1.594 


2.007 


0.5673 


16 


1.589 


2.007 


0.5689 


17 


1.584 


2.006 


0.5703 


18 


1.580 


2.004 


0.5715 


19 


1.576 


2.003 


0.5726 


20 


1.573 


2.001 


0.5736 



series and their accuracy. The discrepancy with the 5-th 
digit of Ref. |23| remains to be explained. 
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